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Abstract 

We investigate the behaviour of the perturbative relation between the photon energy spectrum 
in B — » Xgj and the hadronic P+ spectrum in semileptonic B — > X u iv decay at high orders in 
perturbation theory in the "large-/?o" limit, in which only terms of order a™/?Q _1 are retained. 
The leading renormalon in the weight function W(A, P~) relating the two spectra is confirmed 
to be at u = 1/2, corresponding to nonperturbative corrections at 0(Aqcd /mb)- We show that 
the P-y dependent pieces of the weight function have no infrared renormalons in this limit, and 
so the factorial growth in perturbation theory arises solely from the constant terms. We find no 
numerical enhancement of leading logarithms, suggesting that fixed-order perturbation theory is 
more appropriate than a leadingdog resummation for the extraction of \V u b\. The importance of 
various terms in the expansion of the weight function is studied using a model for the B — > X s ^f 
photon spectrum. Our analysis suggests that higher order perturbative corrections do not introduce 
a significant uncertainty in the extraction of \V u b\. 
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I. INTRODUCTION 



The total rate for the decay B — > X u tv provides a theoretically clean determination of 
the magnitude of the CKM matrix element \y u b\ as a double expansion in powers of a s (mb) 
and AQco/^b PQ- However, to eliminate the background from B — > X c decays, strong cuts 
on the final state phase space are required, which can complicate the theoretical analysis. 
The kinematic regions in which cuts on the charged lepton energy Ee, hadronic invariant 
mass m>x [2] and hadronic light-cone momentum P + = Ex — \ Px I (where Ex and Px are the 
energy and three-momentum of the final state hadrons) [3] are strong enough to eliminate 
the charm background all correspond to the so-called shape function regime, in which the 
local OPE for the partial rate breaks down jH |5] . However, in this region an expansion of 
the partial rate in powers of in terms of non-local operators is still possible, and 

the matrix element of the leading nonlocal operator can be measured in B — > X s 7 decay. 
At leading order in Aqcd/^6, we can write 

dv t = J ^j/M + o^j (i) 

where i labels the decay, Ci(u) is perturbatively calculable, and the shape function f{u) is 
nonperturbative, but universal in inclusive B decays. 1 It is convenient to eliminate the shape 
function altogether, and express integrated rates directly in terms of one another (HI EJ [7] . 
For example, we can write 

f^.f^,A„)^o(y (2) 

where P 1 = m# — 2E^, E 1 is the photon energy and A ~ O(Aqcd)- This defines the weight 
function W(A, P 7 ), which can be calculated in perturbation theory. The 0(AqcdA^&) power 
corrections have been extensively discussed in the literature [El El QUI [HI H2J EES] , and have 
typically been estimated to be below the 10% level for \V u b\ [101 EU EE], although it has been 
argued that sub leading four-quark operators may introduce significant uncertainties [TT] . 

The weight function W(A, P 7 ) has been calculated in fixed-order perturbation theory 
to O(a 2 s (3o) [7]. It is also known to next-to-leading-log order, 0(a" log™ -1 m fe //ij), where 
A*i ~ a/Aqcd^ is the typical invariant mass of the final state [TU] , generalized in [13] . It 

C(u>) can be further factorized into "hard" and "jet" functions; however, for our purposes we will not 
make this decomposition. 
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was shown in [7] that the O(a 2 s (3o) corrections to W(A,P 1 ) are substantial, and the same 
order as the 0(a s ) corrections. Given the size of these corrections, it is important to study 
the convergence of the perturbative expansion. 

In this paper we examine the behaviour of W(A, P 7 ) at higher fixed orders in perturbation 
theory. We work in the framework of the "large-/?o" expansion, in which we calculate all 
terms of order a™/^ -1 [T5l[T6] . While there is no limit of QCD in which these terms formally 
dominate, this class of terms allows us to examine the asymptotic nature of perturbation 
theory, as well as giving an estimate for the size of perturbative corrections. We discuss the 
significance of these terms for the extraction of \V u b\. 



II. BOREL TRANSFORMED SPECTRA AND THE WEIGHT FUNCTION 

Since QCD has an asymptotic perturbative expansion, it is convenient to study the Borel 
transformed series B[R](u) of a quantity R, where 

oo 

R = R-R tICE = J2 r n^ +1 (3) 

n=0 

and 

oo 

*[£](«) =£V. (4) 



n=0 

The expansion for B[R](u) has better convergence properties than the original expansion. 
B[R](u) can in turn be used as a generating function for the coefficients r n 

d n ~ 

r n = —B[R}(u)\ u=0 (5) 

while the original expression R can be recovered from the Borel transform B[R](u) by the 
inverse Borel transform 

POO 

R = R tree + du e~ u/as B[R](u). (6) 
Jo 

Singularities in B[R](u) along the positive real u axis make the inverse Borel transform ill- 
defined. These are referred to as infrared renormalons [T7] , factorially growing contributions 
to the coefficients of the perturbative series, which lead to ambiguities of order (Aqcd/^)™- 
In physical quantities these ambiguities are compensated by corresponding ambiguities in the 
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definition of higher-dimensional nonperturbative matrix elements in the operator product 
expansion of order Aq CD , which render the physical quantity unambigious. 2 

The Borel transform Eq. Q, in the large-/3o limit, may be determined from the order a s 
term, tq, with finite gluon mass following [16J: 

1 /~ . . sin7TM \ 
+ - G (u) ^ . (7) 



Here A is the gluon mass and r ro is a constant. We have used the MS scheme with the 
renormalization scale \x set to the pole mass, m b . The terms Gq(u)/u and arise from the 
renormalization of the graphs involved. 

The weight function W(A, P 7 ) is defined through the relation between the integrated 
B — > Xsj photon energy spectrum and B — ► X u ^z/ charged lepton spectrum, 

r / A x^ f A , p ^ M _ |Kb| 2 vr m| / A ^PWM P^ 77 

u[ ' ~ J + dP + ~ \V tb V^6a em Cf(m b ym b (m b y J ^^^dP, 

+0 ( **°°) (8) 



where A ~ ^qcd in the shape function region, and the normalization is the same as that 
used in 0. Other definitions of W are possible, such as that used in [ID] . As in [7], 
we concentrate on the contribution to the B — > X s 7 spectrum arising from the operator 
7 = (e/16ir 2 )m b s~Lff >J ' 1 ' F^r- While other operators also contribute to the spectrum, for 
the purposes of studying the convergence of the series and estimating the uncertainties from 
higher order terms in perturbation theory we will neglect their contribution and the mixing 
of these with 7 . The factor of m 2 B /ml pulled out of the relation arises naturally, and 
improves the behaviour of perturbation theory for W(A, -P 7 )[Z]. 
Defining the partonic partial rates 

1 dT 77 



T ^ doc 



S(x) +g{x) 



' <R ' ! S(p+) + h(p+) (9) 



T u dp + 



2 Although the renormalon cancellation has only been explicitly shown in some cases in the large-/3o limit, 
it is assumed to hold away from this limit. 
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where T 7 = Gl\V tb V^ 2 a em mf[m b {m b )Cf(m b )} 2 /(32tt 4 ) and T u = G 2 F \V ub \ 2 m 5 b /(192n 3 ) are 
the leading order widths. The partonic variables 

x = 1 — 2E^/m b , p + = (v — q/m b ) • n (10) 

are related to the hadronic variables by 

P 7 = m B - 2P 7 = m b x + A, P+ = E x - \P X \ = m b p + + A (11) 

where A = — m&, g is the momentum of the lepton-neutrino pair, n is a light-like four 
vector in the — q direction and v is the four-velocity of the B meson. Convoluting the 
partonic rate with the shape function to obtain the hadronic rates, we find 

W{A,P y ) = l+ [ A P \h{p)-g(p))dp- [ A Pl g(p)[W(A, P + P,)-l}dp (12) 
Jo Jo 

where the partonic spectra are expanded to leading order in x and p + respectively since in 

the shape function region they are of O^AQCD/mb)- 



Since g(p) and h(p) are 0(a s ) } Eq. (12) may be solved iteratively for W(A, P 7 ). For the 
purposes of this paper, we are only interested in terms of 0(a"/3Q _1 ), for which the last term 



in Eq. (12) does not contribute; therefore, we can write 

W(A, P 7 ) = f£(A - P 7 ) - f? 7 (A - P 7 ) + O«/3 - 2 ) (13) 
where we have defined the integrated partonic rates calculated in perturbation theory, 

f? 7 (A) = 1 [ A d ^dx (14) 



and 

fp u (A) = w T d ^dp + . (15) 



The corresponding quantities W, T^ 7 and are defined by subtracting the tree level con- 
tribution. 

Calculating the parton level photon spectrum with a massive gluon is straightforward, 
and was done in [7j. Integrating the rate with a massive gluon over the endpoint region and 
performing the integral Eq. ([7]), we find the Borel transform of the integrated partonic rate: 
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2 sin nu (1 + u) (3m 2 - 2u - 2) T(u) 2 
nu 2 (m — 2){u — \)uT(2u) cos7TM 

+ - G„(it) . (16) 



It \ 7TM 

Since the operator 7 requires renormalization, the last line arises from the MS counterterm. 
Gq(u) is given by 

oo 

9r. 



= £ ^ n ( 17 ) 

n=0 

and # n are the coefficients of the expansion of Gq(u) [16J 



^ m 2(2m + l)r(4 + 2u) sin™ 
3(m + 2 (m + 1 r(2 + uY nu 

m=0 

The Borel transform of the differential photon spectrum away from the x = endpoint was 
calculated in [T5]. Integrating this result from x = to x = A reproduces the A dependent 



terms of our result, Eq. (16). (The A-independent terms depend on the virtual contribution 
and cannot be directly compared against [18J). 

The calculation of the Borel transform of the semileptonic partial rate T^(A) is signif- 
icantly more involved than for B — > X s 7. The Borel transform of the triple-differential 
B — > X u iv spectrum was calculated in [TH]. Rather than integrate this result over the ap- 
propriate phase space, we instead calculated the integrated rate T U (A) for a massive gluon, 
and then performed the integral ([7]). The result has the comparatively simple form 



8RW1M = e M|fc])(A) 



-2u 



2 5 7 2 1 \ 2sin7ru / A \ u 



3{u — 3) u — 2 3(u — 1) 3m u 2 3(u — 4) J nu \m& / 
T(u) 2 (\ 



+ 7 T7 . — r ( -(9m 4 - 103m 3 - 62m 2 + 38m + 24) 

(m-4)(m-2)(m- l)Mr(2M)cos7TM V3 

-16m(1 + m)(2m- 1) ( , \ + i/j(u) - i/j(2u)))} (19) 

\Sin27TM / 1 J 



where ip(u) = T'(u)/T(u) is the digamma function. 



The Borel transformed weight function is given by the difference between Eq. (19) and 



Eq. (16). Note that the terms proportional to (A/m b ) ju and (A/m&) u sin nu/u , which 



generate the a™ \rt n+l (A/m b ) logs, cancel in the difference. This reflects the universality of 
the leading Sudakov logs. We can resum this contribution by evaluating the inverse Borel 
transform, Eq. ([6]). However the result does not exponentiate because higher powers of logs, 
up to a™ ln 2n double logs, are not included since they are suppressed in (3 . The resummed 
a™ ln n+1 (A/mb) logs from Eq. (Al) and Eq. (A2) are given by 



9i 



_ 4tvu 2 

du e Qs/3 o — 



o Jo 



\m b 



-2u 



A 



Po \ rrib 1 + 2a a s po 1 + 2a 



lit . (1 + a) 
in 



a s (m b )p lvi A 

47T 

The final result for the Borel transformed weight function is 



where a = ^WdEi ln 

47r mb 



(20) 



B[W(A,P 7 )](u) 



T(u) 



2 sin ttu 



(u-5)(3u-4) 



ttu 2 I (u-4)(u-3)(u-2)(u- 1) 



A-P 7 



T(2u)(u - 4)(u - 2)(w - 1) cosvra 



(16(m + 1)(2u-1) (- 

\ V SI 



7T 



sin 27tm 



+ -(5m + 2)(7u + 1) 
o 



1 A ^ / x 2e 5u//3 sin 7tm 
- G («) 



(21) 



where Gq{u) is obtained from Eq. (17) and Eq. (18).Eq. (21 ) is the main result of this paper. 

The Borel transforms can be used to generate the O(a"/?o _1 ) terms in the perturbative 
expansion via the relation Eq. (|5]). Writing 



r^A) = i + ^Q(A)- 



8=1 

oo 



f£(A) = l + ^C^A) 



(4 7r )n 



oo 



W(A,P 7 ) = l + ^H/ n (A,P 7 



i=i 



(4 ?r )n 

' (47r) n 



(22) 



we can easily find the coefficients C*(A), C^(A) and W n (A, P 7 ) to any order. The coefficients 
are given up to n = 5 in Appendix \A\ 



The leading log (LL) and next-to-leading log (NLL) terms in Eq. (A3) are contained 
within the renormalization group resummed NLL result in soft-collinear effective theory 



(SCET), W(A,P 7 )^ T , obtained from j3J QUI [201 121] ■ The SCET result sums lo g s of tf/ m h 
where /if ~ O^KqcD^b)- in the Appendix^ we verify that the leading (3 terms agree with 



Eq. (B2) by expanding in a s {m b ) and then expanding logs of yu?/m| and /xf/(m fe (A — P 7 )). 
Our results also agree with those in [7J |22j [23] . 



III. RESULTS AND DISCUSSION 



A. Renormalons and Borel Resummation 



The leading renormalon ambiguity in both the photon and semileptonic spectra is 
0{A QCD /m b ) due to the pole at u = 1/2 in P[f£(A)](u) and B[f p 77 (A)](u). The diver- 
gence does not cancel between the spectra and gives rise to a pole at u = 1/2 in the Borel 
transformed weight function. This is consistent with the presence of nonperturbative cor- 
rections to W(A, Py) at O^AqcD/n^b) due to subleading shape functions [8]. 

The Borel transform of the weight function can be written in terms of A — P 7 independent 
and dependent pieces, -B[W / ]( M ) an d P[Wi](w) respectively, 

r[ ^i/,a r 5u/3 ( Zsimru f (u - 5)(3u - A) \ 

B[W0](U) ~ 6 {-^{(u-4)(u-3)(u-2)(u-l)- 1 ) 

"pT^y AV T{U t ( ( 16(« + 1)(2« " 1) (-r^- + m ~ V(2«)) 

1 {2u)(u — 4) (u — 2) (u — 1) cos 7tu \ \sm 2iru / 

2, w A\ 1 ( x , s 2e 5M / 3 sinvTM\ . , 

+-(5u + 2) (7u + 1) - - G„(u (23) 

O J J U \ Till J 




B[fc(A,P 7 )](«) = e 5,/3 2sin™ J" " ^ " ^ c I ( ) - 1 1(24) 

L v ' 77U ; ttu 2 (m-4)(u-3)(u-2)(u- 1) 1 1 1 ' 

where we have defined PjWoK^) and B[Wi](it) such that they are finite as u — > 0. Note 
that -B[W / i](m) has no singularities for positive w. Therefore the inverse Borel transform of 



Eq. (24), W±, is well defined and unambiguously resums logarithms of (A — P 7 )/m&. This 
tells us that the poor behavior in the perturbative expansion of the weight function is entirely 
due to the constant terms, Wo, which are generated by P[H / o](w). 

The relevant quantity in determining \V u b\ is the weight function convoluted with the 
B — > X s 7 photon spectrum, as in Eq. (|8|). It is interesting to note that the integrated 



quantity can have a renormalon ambiguity that is not present in the weight function. In 
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order to illustrate this we calculate the Borel transform of Wi, which is renormalon free, 
convoluted with a simple model of the normalized B — > X s 7 spectrum, 

ldT °- * (25) 



r 7 dP 7 r(6)A b ^ 7 

where b = 2.5 and A = 0.77 GeV [10J. This is straightforward to obtain from Eq. (p4j): 



P 



(u) = c W3 2sin7rM (m-5)(3»-4) 



vra 2 (w-4)(«-3)(u-2)(tt-l) 



where r(a, z) = J z °° t a ~ 1 e~ t dt is the incomplete Gamma function. The T(l — u) term in 
Eq. (26) gives rise to a pole at u = 1, which corresponds to an order 0{{KQCD/^b) 2 ) ambi- 
guity in the integrated quantity. This arises because higher order terms in the perturbative 
expansion of W\ have more powers of ln(A — P 7 )/mf, and therefore are more singular near 
the end point. However since the renormalon in PjWo](w) is at u = 1/2, the factorial growth 
in the integrated quantity is dominated by the constant terms in the weight function rather 
than the logarithms. 

It is amusing to notice that if the a™ m(A/mh) n+1 Sudakov logs did not cancel between 
T^(A) and E^A) these would give rise to a ((A — P 1 )/rrii ) )~ 2u term in the Borel transform 




of the weight function. When integrated over P 7 with Eq. (25) this would lead to a pole at 
u — 1/2, the same order as the renormalon in P[W"o](tt). 

Since P[Wi](tt) has no poles in u, the inverse Borel transform of Eq. (24) is well-defined. 



We may therefore use Eq. (24) to sum all terms containing powers of ln((A — P 7 )/mf,) (terms 
of order a™/^ -1 log n ~ m (A — P 7 )/m;,, for n = 1 to infinity and m = to n — 1). While we 
were unable to obtain a closed-form result for this quantity, by expanding Eq. ( |24[ ) in powers 
of u it is straightforward to sum all terms of order «"/3q _1 log n_m ( A — P 7 )/m;,, for n = m + 1 
to infinity and for fixed m > 0, by evaluating the inverse Borel transform 



C F f°° *™ , / / A - P~ 



H/(A,P 7 ) rcsummed = -^J du e'^G m ^u m - x N — — - 1 

CVC m _ lw N ( a s f3 \ m ( ( ^ x a s (3 A - P 7 
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where C m -\ is the coefficient of the u r ' 



■ 1 (A-By\-« 

\ m b J 



term in Eq. (24), and the second line 



follows for m > 0. The constant non-logarithmic terms in the weight function are not 



included in Eq. (27), as they arise from Wo, but may be obtained from the expansion 



Eq. (A3) to give the full resummed logarithmic result. At leading-log (LL), m = 0, we find 



W(A, P 7 )£ = 1 - g In In ^ + 



(28) 



We explicitly show the NLL, the next-to-next-to-leading logarithmic (NNLL) a"/?o 1 In™ 2 
and next-to-next-to-next-to-leading logarithmic (NNNLL) a™/?Q 



?o 1 ln n 3 terms below: 



W(A, P 



NLL 



w(a,pX nll 



a s (m b )CB 

4:71 

a s (m b ) 2 {3 C_ 



14 

3~ 



1 



+ 



(4tt) 2 
/3857 
V 144" 



1 + 6 
/1559 



167 
~36 



2tt^ 
3 



A 216 
16vr 2 



57T 2 

Is" 



[i + by 
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12C(3) 



W(A, P 



NNNLL 
7^/30 



a s (m fe ) 3 ^C f 



(47T) 



/ 65545 
V 3888 



14tt^ 



(i + by 



+ 



/90043 13tt 2 16tt 4 166C(3) 
V 864 ~ 108" ~ ~L5 3 



(29) 



where b 



a s (m b )l3 j n A-P 7 



. These results provide a useful check of our calculation, as they 



47T Oij, 

may be compared with the corresponding resummed expressions in SCET, obtained from 
[31 QUI 1201 EI]- Setting = a/^&(A — P 7 ), we verify that the resummed LL and NLL 



contributions in the large /3q limit, Eq. (28) and Eq. (29), are contained within the RG 
resummed NLL SCET result. 

Finally, the renormalon in the weight function suggests that the dominant contribution 
to its perturbative expansion is from non-logarithmic terms. We can investigate this nu- 
merically by considering the leading logarithmic expansion away from the P 7 — > A end 



point. Following [7], we combine all known terms from Eq. (A3) and Eq. (B2), and take 
the ratio of the various logarithmic terms. While this misses the contributions of terms 
beyond NLL and sub leading in /3q, we can hope that the values below are still indicative of 
the relative contributions of the various terms. Taking mj = 4.8 GeV, a s (m&) = 0.22 and 
/i 2 /m 2 ~ (A — P 7 )/mf, = 1/9 as in [7] we find the following ratios of the logarithmic terms 
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at each order in a,: 



a 3 s : 


0(k>j 


5 3 ) 


0(log 2 ) 


: OQog 1 ) 


: O(log ) 


= 1 : 2.1 


: 1.8 : -6.0 




0(lo { 




0(log 3 ) 


: 0(log 2 ) 


: OQog 1 ) 


: OQog ) 


= 1 : 3.5 : 2.9 : 0.4 : -26 


a 5 s : 


0(lo { 


J 8 ) 


0(log 4 ) 


: OQog 3 ) 


■■ OQog 2 ) 


: ©(log 1 ) 


: OQog ) = 



1 : 4.9 : 4.2 : 1.0 : -2.3 : -119.(30) 

From these results, we can make two observations. First, the renormalon ambiguity in 
the weight function is reflected in the rapid growth of the non-logarithmic terms, which 
dominate the perturbative expansion. However, this bad behaviour of perturbation theory 
is unphysical: in a consistent approach to 0(l/rrib), the renormalon in the weight func- 
tion will cancel with a corresponding ambiguity in the definitions of the subleading shape 
functions. This cancellation would be manifest if the subleading shape functions were con- 
sistently extracted from physical observables, but since they are currently modelled, no such 
cancellation is manifest. We will see in the next section that the estimated uncertainty in 
\V u b\ from the factorially growing terms is small compared to other sources of error, so we 
will not attempt in this paper to absorb the renormalon ambiguity into subleading shape 
functions. These results do, however, underscore the fact that separating the bad behaviour 
of perturbation theory from the 0(l/rrib) corrections is not a well-defined procedure. 



Second, assuming the pattern in Eq. (30) continued to hold beyond the large-/3 and 



NLL terms included here, it indicates that terms which are enhanced by more powers of 
log [if /ml ~ log(l/9) ~ —2 do not dominate over terms with fewer powers of logarithms. 
Since the logarithmic terms do not suffer from renormalon ambiguities, and, therefore, no 
cancellation against the subleading operators is expected for these terms, this pattern should 
not change once subleading operators are consistently included. Thus, these results support 
the conclusion of [7] that fixed-order perturbation theory is more appropriate than a leading- 
log resummation for the extraction of \V u b\ (see also [2U 125]). 

B. Determination of \V u b\ 

From a phenomenological perspective, our results are most useful as an estimate of the 
size of higher-order perturbative corrections to the extraction of \V u b\ via Eq. ([8]). The 
perturbative results in Eq. (A3) for W(A,P 7 ) are plotted in Fig. [I] at different orders in 
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FIG. 1: W(A,i~) from Eq. (A3) is shown including terms up to the following order: (a) O(ag/?o), 



(b) 0{al(3l), (c) O(atPi) and (d) O(a^ 4 ). 



the a"/3o _1 expansion. Throughout this section, we will use the values m b = 4.8 GeV and 
ots{ m b) — 0-22 for numerical evaluations, and take A = m^/m^ = 0.66 GeV, corresponding 
to the kinematic cut which removes the B —>■ X c background. At tree level, the weight 
function is 1 (the dotted line in Fig. [I] and Fig. [2|. Curve (a) in Fig. [I] shows the weight 
function up to O(a;?/?o), calculated previously in [7], while curves (b), (c) and (d) show the 
results to O(c^), 0(a^l) and 0(a 5 s $). 

It is clear from the plots that the perturbative series for W(A, P 7 ) is not converging 
well, as was discussed in the previous section, due largely to the factorial growth of the 
constant terms in W(A, P 7 ). As we will discuss shortly, the results suggest that the optimal 
perturbative estimate is obtained by truncating the series at O(a^), and using the 0(af) 
result as an estimate of the corresponding perturbative uncertainty. In Fig. [2] we therefore 
compare the fixed-order a^/3o result to other perturbative estimates of the weight function. 
Curve (a) shows all known terms up to O(a^): the complete NLL terms from Eq. (B2), 



combined with the additional large (3$ terms in Eq. (A3) that are higher order in the leading 



log expansion. The gray band around the curve gives the perturbative error estimate given 
by the 0{a A s [3l) term. The result is very close to the large-/3o calculation up to O(a^/?o), 
shown in Curve (b). Curve (c) shows the complete NLL resummed result. 

As discussed in the previous section, the integral in Eq. ([8l has a worse perturbative 
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FIG. 2: (a) W(A,P 7 ) with all terms to 0(cq) from Eq. (A3) and Eq. (B2). The grey region is the 



error estimate obtained from the o.^0q term, (b) W(A,P 7 ) up to 0(a^/3g) from Eq. (A3), (c) The 



resummed NLL SCET result, W(A,P 7 )^ T . 

expansion than the weight function itself, since at higher orders in perturbation theory 
W(A, Py) is more singular at the endpoint of integration. Hence, to determine the effects 
of perturbative corrections on the determination of \V u b\, we must look at the perturbative 
expansion of Eq. (|8j rather than that of W(A, P 7 ). For the purposes of estimating the size 
of higher order terms, we adopt the simple model of the normalized B — > X s ^ spectrum, 



Eq. (25). We obtain r u (A), the integrated B — > X u £ui decay rate normalized to the tree 
level value, 

f.(A) = ^ [ A dP + ^ (31) 



u Jo 



shown in Table [Tj We include several more terms than are explicitly shown in Eq. ( |A3[ ) to 
demonstrate that the series appears to converge up to 0(^/3^) and then begins to diverge. 
This suggests that the optimal perturbative result is given by including all terms up to 
O(a^) and using the 0(ot\) contribution to estimate the perturbative uncertainty. At this 
stage, our best estimate of this result is obtained by including all known terms up to 0{o? s ) 



from Eq. (A3) and Eq. (B2), and estimating the uncertainty from the 0{a i s f3 Q v ) term. Table 
III gives r„(A) obtained from the renormalization group resummed LL and NLL weight 



function in SCET, as well as all terms up to 0{a s ) from Eq. (A3) and Eq. (B2). We see 



that the NLL result is in agreement with the optimal perturbative value, within the error. 
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The perturbative uncertainty in \V u b\, estimated from the 0(a A s (3l) terms is approximately 
0.5%, which is far smaller than the order 5% theoretical uncertainty in \V U ^\ from subleading 
shape functions, error in the b quark mass and other sources [TU] . 



IV. CONCLUSIONS 



We have calculated the Borel transform of the B — ► X u £p£ P + spectrum and B — > X s ^ 
P-y spectrum to leading order in Agc^/m?,, from which we determine the Borel transform of 
the weight function. The leading renormalon in W{A 1 P 1 ) is confirmed to be at u = 1/2, 
corresponding to nonperturbative corrections at 0(AQcr>/mb). The a™/?^ -1 terms are easily 
obtained from the Borel transform of the weight function and are given analytically to 
n = 5. We are able to resum logarithms of (A — P 7 )/m;, in the large (3q limit of the weight 
function since the relevant terms in B[W(A, P 7 )](u) are renormalon free. However we show 
that integrating these terms over P 7 introduces a renormalon. Comparing all known terms 
in the perturbative expansion of the weight function, we find no numerical enhancement 
of leading logarithms, suggesting that fixed-order perturbation theory is more appropriate 
than a leading-log resummation. 

From our results we estimate the size of higher-order perturbative corrections on the 
extraction of \V u b\ using a model for the B — > X s 7 photon spectrum. We have shown that 
T M (A) begins to diverge beyond 0{a A (3l) in the (3q expansion of the weight function. This 
suggests that the best perturbative estimate is given by including terms up to O(a^) with 
the theoretical uncertainty given by the a^/So term. We show that this result is in good 
agreement with the resummed NLL SCET result. 
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Tree 


0(a s ) 


o(« s 2 A)) 


O(«?/3 2 ) 


0(ajf3i) 




0(a«$j) 




1 


1.08 


1.15 


1.17 


1.16 


1.12 


1.04 


0.88 



TABLE I: T U {A) for different orders in the "large-/3 " expansion of W(A,P~), Eq. (A3) 
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Tree 


SCET LL 


SCET NLL 


All known terms to 0(a 3 ) 


1 


1.10 


1.18 


1.17 



TABLE II: T U (A) for the resummed LL and NLL weight function in SCET, and all terms up to 
0(a 3 s ) from Eq. ^ and Eq. rigg}. 
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APPENDIX A: EXPANDING OUT THE FUNCTIONS 



A A An 2 

CHA) = -21n 2 — -71n— - — -5 

m b m b 3 

™/an n, 3 A 13. 2 A / 85 2tt 2 \ A 91vr 2 631 

C|(A = 21n 3 — + -ln 2 — + -- + — In 4(3- — - — 

m b 6 m& \ 6 3 / m b 18 36 

_, /AN 7 4 A 1 ,A /275 2tt 2 \ i2 A 1, 2N , A 

CJ(A =- In 4 — + -ln 3 — + — -— In 2 — + - -581 + 587T 2 In — 
3 m fc 3 m b \ 18 3 / m b 18 m b 

+ 3^4 (~ 12727 ~ 6366?r2 - 1087r4 - 13824C(3)) 

ns/M 01 5 A 35 4 A / 35 2tt 2 \ 3 A /6029 29tt 2 \ 2 A 

CM A) = 3 In 5 In 4 h 1 In' 3 h In 2 — 

4V ' m b 12 m b V 2 3 / m b \ 108 6 J m h 

( 9557 235vr 2 vr 4 \ , A A/ , 555 A/ s ,/ 24959 

+ + In 72C 5 C 3 + 7T 2 8C 3 

V 108 18 5 J m b sv ; 2 su \ 324 sv ; 

57tt 4 283555 



10 2592 

^ S/AN 62 6 A 33 5 A /395 2vr 2 \ 4 A / 2543 58vr 2 \ 3 A 

CI(A) = In 6 1 In 5 h In 4 M 1 In 3 — 

5V ; 15 m b 5 m 6 V 18 3 / m 6 \ 27 9 ) m b 

/32171 235vr 2 2tt 4 \ , 2 A / 50189 4429tt 2 29vr 4 \ , A 



162 81 15 / m b 

7392583 154997vr 2 3932tt 4 494tt 6 7452 
19440 486 75 315 5~~ ^ ' 

/ 205219 496tt 2 >/ A >/ . /A . 
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m b o 

, A o A 

2 In 3 — + 3 In 2 — 

m b m b 



3 mf, 
168313 



m b 
- + 
A 



36 



113 2tt 2 



- - ln° — + 
9 m b 

2135tt 2 7tt 4 



359 
l8~ 



27T 2 

IT 



taA_ 16C(3 )-i^ + l^ 
m 6 SV ; 6 144 

5045 34tt 2 



In 2 — + 



108 



9 



In 



m b 



- 98C(3) 



2592 108 5 

3 i n s A _ 5 ln 4 A + I (_i 33 + 4^) ln 3 A + 

m b 2 m b 6 m b 



+ 



180229 319tt 2 
+ 



1296 



18 



7T 

~5 



In 



m b 



432C(5) 



/ 16735 17tt s 
V~216 3~ 
1807 



In 2 



m b 



+ 40vr 2 C(3) 



13129tt 2 
432 



79tt 4 11428313 
+ 



6 



15 m?, 
/21581 



31104 
A 94 



+ 



+ 



V 72 
/ 29741 

V 54 



15 m& 
319tt 2 2tt 4 

1408tt 2 



479 2tt 2 
~18~ ~ ~3~ 

In 2 h 



m b 



2215 68vr' 



1? 



668117 13535tt 2 



- 672 C (3)j C(3) + 



1296 
8231tt 2 



48 



162 
2774tt 6 
315 



9 

34tt 4 
~L5~ 
1649tt 4 
30 



ln d 



A 

m b 
A 



+ 



In 4920C(5) 

m b 

64526377 



31104 



(A2) 
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IVl (A, E() = __ ln __7__ + _ 

^(A, P 7 ) = | .„> ^ - H ln ^ + ^ _ «!£ _ 12C(3) 

6 nib o nib 144 9 

„, , A „ , 5 3 A - P 7 14 2 A - P 7 /5tt 2 1559 \ , A - P 7 90043 

W 3 (A, P 7 ) = — hr 1 H In 2 1 + ln 1 H 

v ' 7; 9 m b 3 m b \ 9 108 / m b 864 

13tt 2 16vr 4 166 



C(3) 

5 , 4 A-P 7 14, ,A-P, /1559 57r 2 \ 2 A-P 



108 15 3 
jy 4 (A,P 7 ) = --In 4 — ^-^ln 3 ^^+ ( ) ln 



/1559 5vr 2 \ 
\~T2 6~) 



12 mj 3 nib \ 72 6 J rrib 

14vr 2 65545 \ , A - P 7 . . /71 A w s 

ln * - 360C(5) - — + 32vr 2 C(3) 



112vr 4 60449tt 2 14830973 
+ + 



15 1296 31104 

rir , . „ , 1 , 5 A - P 7 14 , 4 A - P 7 / 1559 10tt 2 \ , 3 A - P 7 
W 5 {A, P 7 ) = -- In 5 1 + — In 4 ^ + — + — — In 3 



3 nib 3 nib \ 54 9 / m b 

/ 65545 28tt 2 \ „ A - P 7 / 266605 1559vr 2 tt 4 \ A - P 7 



/ 65545 28tt 2 \ 2 A - P 7 / 266605 1559tt 2 tt 4 \ 
V~648 3/ 11 m fe + V 1296 + ^4 ~3~y 11 

lZli5 C(5 ) + r_576C(3) + ^ _ 3 04^ C (3) - ^ ^ 



10 J " 7 21 50 

1906687tt 2 381772549 

H 1 . (A3) 

3888 155520 v ' 



APPENDIX B: THE WEIGHT FUNCTION TO NLL ORDER 

The renormalization group resummed NLL weight function has been calculated in SCET, 
PES El [21]. By expanding W(A, P 7 )^ T in a s (nib) and re-expanding the logarithms of 
/i 2 /m 2 and /x 2 /(mb(A — P 7 )) we find 

W(A, P 7 & = i + J2 W n (A, P 7 )^ T a ' ( y F (Bl) 



i=i 
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and the first coefficients are given by: 
^(A,P 7 ) SCET - -gin— - + — 



W 2(A , P 7 )- = (*J* + 92\ A-R, /_ 14ft 128 85^ _ 5122 \ A-R, 

2V ' 7;scET V 6 27 J m b \ 3 3 sv ; 27 81 J m b 



W 3 (A, P 7 ) St = ( - 5 4 - ^ - ^ ) in' + f M + C _« C(3 ) - & ' 



9 27 243 J m b \ 3 V 3 27 

11501 \ „ 2560 512vr 4 2392tt 2 68155V 2 A -P~ 
+ -TT^ A - -R=-C(3) " — 7i T^T" + "^77" ln 



162 /' u 27 ^ w 135 243 162 / m 6 



W(A p \nll ^5/? 3 253/3 2 808ft 27584 \ 4 A - P 7 / 14/? 3 f 128 



175tt 2 19981 \ ^ 2 /2560 A/ON 512vr 4 3220tt 2 243991 \ . 

^ oTo^ $ + ~^C(3) + -^r + j^r- A 



51 243 



65536 47104 w „, 2048tt 4 56560tt 2 21384356 V 3 A - P 7 



+ C(5) + C(3) + + ln d 

27 si ) 243 si ) 243 21g7 6561 J 



m 6 



,nll f Po 230^0 2828$ 55168/3 462080 V 5 A - P 



W 5 (A,P 7 )^ T = - - - - In 



7 



3 81 243 2187 19683 / m b 

jWS.f 32 557T 2 , 89585 \ 3 / 7040 1408tt 4 

+ V 3 + V 3 adj 27 + 972 J ^ + { 81 405 
111327T 2 548459\ n2 ( 32768 >/ , 23552 ^ , 1024vr 4 

-^r + -iwr J * + (-^ c(5) - ~iV (3) - 

35552tt 2 3675094 \ n 1310720V x 827392 >/ , 262144vr 6 
+ j ft 77^C(5) - ^T^C(3) 



729 729 V 243 2187 15309 

188416tt 4 1213696tt 2 1371073480 \ , 4 A - P 7 

1 ln -. (B2) 

10935 19683 59049 J m b K ' 



We verify that the leading ft terms agree with Eq. ( A3 ) . 
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